BOHMIAN MEASURES AND THEIR CLASSICAL LIMIT 
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Abstract. We consider a class of phase space measures, which naturally arise 
' in the Bohmian interpretation of quantum mechanics. We study the classical 

limit of these so-called Bohmian measures, in dependence on the scale of os- 
(-H , cillations and concentrations of the sequence of wave functions under consid- 

eration. The obtained results are consequently compared to those derived via 
|__^ ' semi-classical Wigner measures. To this end, we shall also give a connection 

to the theory of Young measures and prove several new results on Wigner 
measures themselves. Our analysis gives new insight on oscillation and con- 
centration effects in the semi-classial regime. 
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Prob x(t)60 = / H £ (t,x)\ z dx, 



1. Introduction 

In this work we consider a quantum systems described by a time-dependent 
wave-function i/j s (t, •) G L 2 (M. d ; C) . The dynamics of -0 e is governed by the linear 
Schrodinger equation 

(1.1) ied t 4> e = -yA^ e + F(x)V> e , V> £ (0,x) = ^(x). 

where x G R d , (£M and V = V(i) a given real-valued potential. Here, we already 
rescaled all physical parameters, such that only one dimensionless parameter e > 
remains. We shall mainly be interested in the classical limit e — > 0+, and so in the 
following refer to e as the semi- classical parameter. In the usual interpretation of 
quantum mechanics, the wave-function yields a probabilistic description of the 
position of the particle X(i) G R d , at time t G R. More precisely, 

m 

is the probability of finding the particle a time t G K within the region f2 C R d . 
This requires the wave function to be normalized ||V' e (^j OIIl 2 = 1- 

In more generality, we recall that, although the wave function ip e itself is not a 
physical observable, (real-valued) quadratic quantities of ip e yield probability den- 
sities for the respective physical observables. Two important examples of such den- 
sities, describing the expected values of observables (in a statistical interpretation), 
are the position and the current density, i.e. 

(1.2) p £ {t,x) = \ip £ (t,x)\ 2 , J e (t,x) = elm(iF(t,x)V^ £ {t,x)). 

It is easily seen that if ip £ solves (jl.ll) , then the following conservation law holds 

(1.3) d t p e + div J e = 0. 

Similarly, one can define the total energy of the particle, which is conserved along 
sufficiently smooth solutions to (|1.1[) . In our case it is given by 

(1.4) EW(t)\ = ^-\ \Vr(t,x)\ 2 dx+ [ V(x)\r(t,x)\ 2 dx, 

1 JR d JR d 

i.e. by the sum of the kinetic and the potential energy. 

The semi-classical regime of quantum mechanics corresponds to situations where 
e«l. Note that e corresponds to the typical wave-length of oscillations within the 
sequence of wave functions i\) e . In view of (jl.l[) . this is a highly singular asymptotic 
regime and thus analyzing the limiting behavior of expectation values of physical 
observables requires analytical care. In particular, the limit of the highly oscillatory 
wave function ip E itself is of almost no relevance due to the non-commutativity of 
weak limits and nonlinear functions. 

The conservation law (|1.3[) . is also a possible starting point of the Bohmian 
interpretation of quantum mechanics [9] (see also j!4j for a broader introduction). 
To this end one introduces the velocity field 

(1-5 u £ (t,x) := =glm 

p e {t,x) V ^ (M) 

which is well-defined, expect at nodes, i.e. zeros, of the wave function tp E . Ignoring 
this problem for the moment, the Bohmian dynamics of quantum particles is gov- 
erned by the following system of ordinary differential equations for the macroscopic 
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position vector: 

(1.6) X E (t, x) = u e {t, X e {t, x)), X e (0, x) = xe R d . 

This can be considered as the Eulerian viewpoint of Bohmian mechanics, with 
u e the associated Eulerian velocity. The solution on the Schrodinger equation, 
usually called pilot-wave, is hereby used to define a dynamical system governing 
the time-evolution of particles via (jl.6l) . In addition one assumes that initially 
the particle's position is not completely known, but described by the probability 
distribution p e (0,x) = |^>g(a;)| 2 . This probabilistic feature of Bohmian mechanics 
can be understood as a lack of knowledge about the fine details of the considered 
quantum mechanical system, analogously to the situation in classical statistical 
mechanics, cf. |14) . 

In order to study the semi-classical limit of Bohmian mechanics, one might hope 
to analyze directly the limit of the trajectories X £ (t) as e —> 0+. To our knowledge, 
the only rigorous result in this direction has recently been given in |15j , where the 
authors restrict themselves to the case of so-called semi-classsical wave packets, see 
Section 5.3. In the present work we shall consider a more general situation which, 
as e — > 0+, leads to concentration and oscillation effects (within the particle and 
current density) not present in the case of |15j . In particular it will become clear 
from the Examples given in Section 5 and the comparison to WKB asymptotics 
(see Section 6) that a direct study of the limiting Bohmian trajectories is very hard 
and a complete solution of this problem seems out of reach so far (cf. [5] where 
the most recent results on the limiting behavior of singularly perturbed dynamical 
systems are proved, albeit in a setting much simpler than ours). Instead, we shall 
pass to the corresponding Lagrangian point of view on Bohmian dynamics and 
argue that this naturally leads to the introduction of a certain class of probability 
measures on phase space, which we shall call Bohmian measures. These measures 
are supported on sub- manifolds (of phase space) induced by the graph of the velocity 
field u E (t, x). They consequently evolve via the Bohmian phase space flow and can 
be shown to be equivariant with respect to this flow. In addition the first and 
second moment of these measures yield the correct quantum mechanical position 
and current density. It is therefore natural to consider the classical limit of Bohmian 
measures in order to analyze the emergence of classical dynamics from Bohmian 
mechanics (see also Remark 13.81 for some preliminary observations regarding the 
connection to the classical limit of Bohmian trajactories). 

The main analytical tool for studying the classical limit of Bohmian measures 
will be the theory of Young measures, which provides information on weak limits of 
oscillatory sequences of functions. The obtained limit will then be compared to the 
well established theory of (semi-classical) Wigner measures associated to ifj £ (t), see 
e.g. [24l [19j [31] . These are phase space measures which, after taking appropriate 
moments, are known to give the correct classical limit of (the probability densities 
corresponding to) physical observables. We shall prove that the limiting Bohmian 
measure of ip £ (t) coincides with its naturally associated Wigner measure locally 
in-time. That is, before caustic onset, where the first singularity occurs in the 
solution of the corresponding classical Hamilton- Jacobi equation, see Section 6. 
Furthermore, we shall argue (by examples) that in general the limiting Bohmian 
measure differs from the Wigner measure after caustic onset. In the course of this 
we shall also prove new results on when Wigner transforms tend to mono-kinetic 
Wigner measures. 
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The purpose of the present work is thus twofold: First, to gain some information 
on the classical limit of Bohmian mechanics (see in particular Section 6). Second, 
to give further insight on oscillation and concentration effects in the semi-classical 
regime by means of two physically natural, yet again mathematical very different, 
descriptions via phase space measures. Moreover, we believe that our analysis 
may very well be used as a first building block towards an optimal transportation 
formulation of quantum mechanics, by combining our results with those given in 
HH7] and [3D], see also Remark below. 

2. A Lagrangian reformulation of Bohmian mechanics 

2.1. Existence of Bohmian trajectories. We start with some basic assumptions 
on the potential V. Since in this work we shall not be concerned with regularity 
issues we assume 

(A.l) yeC°°(M d ;M), V(x)^0. 

This is (by far) sufficient to ensure that the Hamiltonian operator 

(2.1) H* = -jA + V(x), 

is essentially self-adjoint on D(H £ ) = C%° C L 2 (R d ;C), cf. [3D]. Its unique self- 
adjoint extension (to be denoted by the same symbol) therefore generates a unitary 
Co— group U E (t) — e~ ltHe l E on L 2 (R d ), which ensures the global existence of a 
unique solution ip £ (t) — U £ (t)ipo of the Schrodinger equation (jl.lj) . such that 

\\ P e (t,-)\\ L , = \\r(t,-)\\ 2 L2 = \\r \\h, VieM. 

From now on, we shall also impose the following assumption on the initial data: 
(A.2) Vo G C°°(R d ), with H^olU 2 = 1 and E Wo\ < uniformly in e. 
Since U E (t) and H e commute, we also have that the total energy is conserved, i.e. 

#(t)]=S[4 VteR, 

and thus, in view of (|A.2p . E[ip e (t)] is uniformly bounded as e — » + and for all 
(el. In addition, the dispersive properties of U e (t) together with the assumption 
dH]| imply that if Vo G C Q (R d ), for a ^ 0, then V £ (V) G C a (W l ) for all times 
t € R. In the following, we denote 

H/liri :=\\f e \\^ + \\eVf s \\ L 2 

and we say that a sequence f e = {/ £ }o<e^i is uniformly bounded (as e — > 0+) in 
H^(R d ;C), if 

sup ||/ e ||m < +oo. 

0<£^l 

Note that the two conservation laws given above, together with (jA.ip , (|A.2[) , imply 
that ip £ (ti ') G H^(M. d ) uniformly bounded as e — >• 0+ and for all t € R. Moreover, 
in view of (|1.2p we have 

(2-2) \\r(t,-)\\ L1 ^e\\Vi/(t,-)\\ L 2\\^(t,-)\\ L 2^E[r Q }, 

and we conclude that for all t £ R: J £ (t, ■) e L 1 (R d ;R d ) uniformly as e — > O-i- , 
provided assumptions (jA.ip and (|A.2I) hold. 

Next, we recall the main result of [35] (see also [7]) on the global existence of 
Bohmian trajectories. 
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Proposition 2.1. Let (|A.1|) , (lA.2j) be satisfied. Then the map Xf : x H> X s (t, x) £ 

M. d induced by (11. 6[) exists globally in-time for almost all x £ M. d , relative to the 
measure p$ = \ipQ(x)\ 2 dx and Xf £ C 1 on its maximal open domain. 

Moreover, the probability density p £ (t, •) is the push-forward of the initial density 
Pq under the map X\, i.e. 

p s (t)=X £ #p £ . 

We consequently infer that the Bohmian trajectories, defined through the ordi- 
nary differential equation (|1.6[) . exist pg — a.e. and that for any compactly supported 
test function a £ C (IR d ) it holds 

(2.3) / a{x)p £ {t,x)dx = a{X £ {t, x))p £ Q (x)dx. 

This property is also called equivariance of the measure p £ (t, •) in 0132]. In addition 
we may interpret (|2.3|) as a way of giving sense to the solution of the continuity 
equation 

(2.4) d t p £ + div(p £ u £ ) = 0, 

where u e is given by (|1.5[) . Due to the possible occurrence of nodes in ip £ (t,x), the 
vector field u £ {t,x) is in general not Lipschitz in x. In fact, not even the general 
existence theory [3j [13] for velocity vector fields which only have a certain Sobolev 
or BV regularity applies to Bohmian trajectories. The property ()2.3[) therefore can 
only interpreted as a very weak notion of solving the continuity equation (|2.4[) . 



2.2. Bohmian measures on phase space. We shall now reformulate Bohmian 
mechanics, in its (well-known) Lagrangian formulation. To this end we first intro- 
duce the Lagrangian velocity 

P e (t,x) = X £ (t,x), 

for which we want to derive an equation of motion. In view of (jl.6p . we can 
differentiate P £ (t,x) p e — a.e. to obtain 

P £ {t, x) = d t u e {t, X £ (t, x)) + (X £ (t, x) ■ V)u £ (t, X £ (t, x)) 

= d t u £ {t,X £ (t,x)) + (u £ {t,X(t,x)) ■ \7)u £ {t,X £ (t,x)). 

To proceed further, we need an equation for the velocity field u £ . To this end, we 
recall the well-known hydrodynamic reformulation of quantum mechanics, where one 
derives from (|1.1[) a closed system of equations for the densities p £ , J £ . Assuming 
that ip £ is sufficiently differentiable, they are found to be (see e.g. [TBI ) 

d t p £ + div J £ = 0, 

(2.6) { fJ £ J £ \ e 2 f&JW 



p £ J 2 \ V? 

Under the regularity assumptions on V and ip £ stated above, the weak formulation 
of the quantum hydrodynamic equations (|2.6I) holds in a rigorous way, i.e. each of 
the nonlinear terms can be interpreted in the sense of distributions, see [16j Lemma 
2.1]. 

Remark 2.2. Let us also point out that the hydrodynamic picture of quantum 
mechanics originates in the seminal work of E. Madelung [25) . who interpreted 
p £ , J £ as a description of a continuum fluid instead of a single particle. 
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Identifying the current as J £ — p £ u £ we can formally derive from (|2.6j) the 
following equation for u e : 

(2.7) d t u £ + (u £ • V)u e + W = %- V ^ A%/ ^ 



2 V 

The right hand side can be seen as the gradient of the so-called Bohm potential (or 
quantum potential), given by 

e 2 



(2.8) V £ B 



2 sf(F 

Plugging (|2.7|) into (|2.5p we finally arrive at the following system of ordinary dif- 
ferential equations: 

' X £ = P £ , 



(2.9) 
subj( 

the initial velocity given by 



P £ = -W(X £ ) - VV§(t,X £ ), 
subject to the initial data X £ (0,x) = x and P £ (0,x) = u £ (0,x), where u £ (0,x) is 

u £ {x) = elm » 

Remark 2.3. Note that the system (|2.9p fully determines the quantum mechan- 
ical dynamics. It can be regarded as a system of ordinary differential equations, 
parametrized by the spatial variable x g R d through the initial data, where the 
position density p £ (and its derivatives up to order three) have to be determined 
additionally. At least numerically, p £ can be computed via ray tracing methods or 
particle methods [121 GSl E5] , based on the push- forward formula (|2.3|) . i.e. evolving 

p £ (t, x) « ^2 c « s ( x - x£ ( t ^ 

Strictly speaking, though, this requires the trajectories X £ (t, x) for all x £ R d 
at time f e I. However, it approximately yields the solution of the continuity 
equation (|2.4I) . without solving the Schrodinger equation. This fact makes the 
Lagrangian reformulation interesting for numerical simulations, in particular in 
quantum chemistry, see e.g. [2TJ [26l [34] for applications and [35] for a general 
overview. However, we caution that the Bohm potential Vb is singular at nodes 
of the wave function, which generates significant numerical difficulties in actual 
computations (cf. [12 for more details). 

In order to give (|2.9[) a precise mathematical meaning we shall in the following 
introduce what we call Bohmian measures on phase space R^ x R^. To this end we 
denote by A4 + (M. d xRj) the set of non- negative Borel measures on phase-space and 
by (•, •) the corresponding duality bracket between A4(R^ x R^) and Co(R d x R^), 
where Co(R^ x R^) is the closure (with respect to the uniform norm) of the set of 
continuous functions with compact support. 

Definition 2.4. Let e > be a given scale and ip £ G H^(M. d ) be a sequence of 
wave functions with corresponding densities p £ , J £ . Then, the associated Bohmian 
measure j3 £ = /3 £ [ip £ ] G M + {R d xljj) is given by 

{P £ ,V)--=I p E (x)cp(x,^l)dx, V^Co(R'x<). 
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Note that in the definition of j3 £ a fixed scale e is imposed via the scaling of the 
gradient in the definition of the current density (|1.2j) . Formally, we shall denote 
the Bohmian measure by 



(2.10) F( x ,p)= p *( x )5[p--^ j ee W(x)\'S (p- slm v ( . ( (j 

where S is the usual delta distribution on M. d . Obviously, (|2.10l) defines a continuous 
non-negative distribution on phase space. In addition, the first two moments of /3 £ 
satisfy 

(3 e (x, dp) = p £ {x), [ p p £ (x, dp) = p e {x)u e (x) = J s {x). 

However, higher order moments of /3 s in general do not correspond quantum me- 
chanical probability densities (defined via quadratic expressions of ip s )- I n partic- 
ular, the second moment of /3 e yields 

\n\ 2 1 
^-0 B (x,dp) = -p e (x)\u e (x)\ 2 . 



In classical kinetic theory, this would be interpreted as the kinetic energy density 
of the particle. However, in view of 

(2.11) E kin m := £ 4 [ Wf(x)\ 2 dx=U HM! + ^/ \V^W\ 2 dx, 
2 Jri 2 J Rd p £ (x) 2 J Rd 

we see that the second moment of (3 e is not what in quantum mechanics would be 
called a kinetic energy density, since it does not account for the second term cx e 2 . 
Note that this term formally goes to zero in the classical limit e — > 0+. 

To proceed further, we shall introduce the following mapping on phase space, 

(2-12) *f : (x,p) ^ (X^t,x,p),P s (t,x,p)) 

where X s , P e formally solve the ODE system (|2.9I) for general initial data x,p € M. d . 
Note that this <J>f is not necessarily well defined as a mapping on the whole phase 
space. However, from Proposition ^. 1[ it is straightforward to conclude the following 
existence result. 



Lemma 2.5. Under the same assumptions as in Proposition \2.1\ the mapping $j 
exists globally in-time for almost all (x,p) £ M. 2d , relative to the measure 

Po(x,p) = Po( x ) S (P - u o( x ))- 
Moreover $j is continuous w.r.t. ig I on its maximal open domain and 

(3 £ (t)^^ t #(3 E . 

Proof. First note that <&f when restricted to {graph(ug)} C R d x R d is well defined 
/3q — a.e., since the map Xf established in Proposition 12 . 1 1 does not run into nodes 
of ip £ (t, ■) for almost all x relative to pg. Now, let us w.r.o.g. consider test function 
tp(x,p) = a(x)x(p) € C b (M. 2d ) and denote u e = J e / p e . Then, we have 



(/3 E (t),<p} = I a(x) X (u e (t,x))p £ (t,x)dx 

<i(X e (t, x)) X (u E (t, X E (t, x))) Po (x)dx, 
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where for the second equality we have used (|2.3[) . By definition u £ (t, X £ (t, x)) — 
P e (t,x), hence 



Since P £ (0, x) = u £ (0, X £ (0, x)) — Uq(0, x) the assertion of the lemma is proved. □ 

Thus (3 £ {t) is supported on a subset of phase space given by the graph of the 
velocity field and evolves through the quantum mechanical trajectories, induced 
by (I2.9|) . in the sense of a push- forward for measures. This result is conceptually 
important, since it established the existence of a i-parametrized family of phase 
space measures fi £ {t) = /3 £ (t,x,p), which can readily be compared to the concept 
of Wigner functions w £ (t,x,p), see Sectional and which encodes the full quantum 
mechanical dynamics of p £ (t, x) and J £ (t, x). It is therefore natural to consider the 
limit of (3 £ (t) as s — > 0+ in order to gain insight into the classical limit of Bohmian 
mechanics. This will be the main task of the upcoming sections. 

Remark 2.6. Lemma [2.51 formally allows to interpret f3 £ (t) as a solution of the 
following nonlinear kinetic equation 



subject to initial data /3q as given in Lemma |2~51 As before, equation (|2.13[) can be 
seen as a conservation law in phase space (endowed with a complex structure) . The 
main problem of (|2 . 1 3[) is, that the term W x Vg ■ V p /3 e can not be defined in the 
sense of distributions in a straightforward way. We remark, however, that in the 
purely diffusive setting of so-called quantum drift diffusion models, this mathemat- 
ical difficulty was overcome by using Wasserstein gradient-flow techniques [20] . We 
believe that a combination of [3U] with the results given in [5J [T7] can lead to rig- 
orous mathematical results on (|2 . 1 3[) . opening the door for a new interpretation of 
Bohmian mechanics via optimal mass-transportation. From a mathematical point 
of view, the study of (|2.13j) is also interesting for general measure valued initial 
data, even though, the connection with the Schrodinger equation is lost in such a 
situation. 



We recall that the assumptions on the initial wave function ipQ together with the 
arguments given at the beginning of Section [2] imply that for all t € M the solution 
of the Schrodinger equation ip £ (t) is uniformly bounded in H^(R d ) as e — > + , with 
a bound independent of time (namely, the initial energy). Since the latter will 
be the main technical assumption needed from now on, we shall for the sake of 
notation suppress any time-dependence in the following and formulate results on 
Bohmian and Wigner measures associated to general sequences of L 2 functions ip 
with uniformly (in e) bounded mass and energy. In Section 6 we shall get back to 
time-dependent setting and connect it to the results of the previous sections. 




a(X £ (t, x)) X {P £ (t, x))p (x)dx = <$* # p £ S p=P e (0) , ax) ■ 



(2.13) d t (3 £ + p ■ V x (3 £ - V X (V + V§) ■ V p /3 £ = 0, 




3. The classical limit of Bohmian measures 



3.1. Existence of limiting measures. We start with the following basic lemma, 
which ensures existence of a classical limit of f3 £ . 
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Lemma 3.1. Let ip £ be uniformly bounded in L 2 (R d ). Then, up to extraction of 
sub-sequences, there exists a limiting measure (3° = (3 6 A4 + (R d x R d ), such that 

(3 e £ ^ (3 in M + (Ri xR d )w- *. 
Proof. In view of Definition 12.41 we have, for all test functions ip S Co(R d xlj): 

l(/3 E ,^)l < II^IU»(E 2 «»)l|p e ||Li(K<i) < +00, 

uniformly in e, by assumption. By compactness, we conclude that there exists 
a sub-sequence {£„}„£N, tending to zero as n — > oo, such that (3 En " — V (3 in 
M + {Ri x R d ) weak - *. □ 

When ip e — ?p e (t) evolves according to the Schrodinger equation with initial 
data such that ||V>oIIl 2 = 1, then f3 e is in L°°{R t ,M + (Ri x Rf)) uniformly as 
e —> 0+, since p e is in L oa (R t ,L 1 (R d )) uniformly as e — > 0+. Thus there exists 
a sub-sequence (which we shall denote by the same symbol) and a t-parametrized 
family of limiting probability measures f3° — f3°(t), such that (3 e tends to f3° in 
L°°(R t , M + (R d xlp)) weak - *. This is of importance when we shall get back to 
Schrodinger wave functions in Section 6, in particular for Proposition 6.1. 

Next, we shall be concerned with the classical limits of the densities p e , J £ . Since 
they are both uniformly bounded in L 1 (R d ), provided ip £ is uniformly bounded in 
H^(R d ), we conclude, that, up to extraction of a subsequence, it holds 

(3.1) p e £ 4 + p, inAi + (^;t)w-*, J e e =^ + J, in M + (R d x ; R d ) w - *. 

Moreover, it has been proved in [16], that J <C p in the sense of measures and thus, 
by the Radon-Nikodym theorem there exists a measurable function u, such that 

(3.2) dJ = u(x)dp. 

Formally, the function u(x) G R d can be interpreted as the classical limit of the 
Bohmian velocity field u £ . The following statement gives the connection between 
the limits (p, J) and f3. 

Lemma 3.2. Let i\) e be uniformly bounded in H^(R d ). Then 

(3.3) p{x)= [ 0(x,dp), J(x)= ( P P(x,d P ). 
Moreover, we also have 

(3-4) lim ff (3 £ m(dx,d P )= [[ P(dx,dp), 

provided that the sequence ip e is compact at infinity, i.e. 

lim lim / \ip £ (x)\ 2 dx = 0. 

Thus, the classical limit of the densities p e , J e can be obtained from the limiting 
Bohmian (phase space) measure (3 by taking the zeroth and first moment. In 
addition no mass is lost during the limiting process at |x| + \p\ = +00, if in addition 
ip £ is compact at infinity. 

Remark 3.3. Note that the property of ip £ being compact at infinity [HO HI] can 
be rephrased as \ip e \ 2 being tight. 



10 



P. MARKOWICH, T. PAUL, AND C. SPARBER 



Proof. We first prove assertion (|3.3|) for p £ . To this end let a £ Co(R d ) and write 
(j(x)(3 e (dx,dp) = If a(x)xR.(p)l3 e (dx : dp) 

a(x)(l - XR{p))P e {dx,dp), 



where for a given cut-off R > 0, xr S Co(K d ), such that: ^ XR ^ 1 an d Xr(p) = 1 
for < i?, as well as xip) = f° r \p\ > R + 1. In view of Lemma [3.11 the first 
integral on the r.h.s. converges 

o-{x)XR,(p)P e (dx,dp) // a(x)xR{p)(3(dx,dp). 

■ J jR2d 

On the other hand, the second integral on the r.h.s. is 

v( x )i l ~ XR{p))l3 £ {dx,dp) < / ^ e (£c)o-(i)l{| tt s| >Ji }da; ) 

R 2d JR d 

where u e = ^ and In denotes the indicator function of a given set ft C R d . We 
can now estimate 

If C 

p e (x)a(x)l { \ u e\ >R} dx < — / |J £ (i,x)|cr(x)dx ^ — , 

tt J R d H 

where C £ R+ is independent of e. Here, the last inequality follows from (|2.2j) 
together with the uniform bound of i\) e in H^(R d ). We can therefore take the 
respective limits e — ¥ 0+ and i? — > +oo, to obtain the desired statement for the 
position density p £ . The assertion (|3.3I) for J e can be shown analogously. Finally, 
in order to prove (|3.4p , we refer to [T9l [24] , where it is shown that 

lim / p £ (x)dx— I p(x)dx, 
£ -* + Js, d J& d 

provided p e = \ip £ \ 2 is tight. Jointly with (|3.3[) . this directly implies (|3 .4[) . □ 

The main task of this work is henceforth to study the limit /3 S Ai + (R d . xRj). 
In particular, we want to understand under which circumstances (3 is mono-kinetic. 

Definition 3.4. We say that (3 £ A4 + (M. d , x M~) is mono-kinetic, if there exists a 
measure p € A / 1 + (M^) and a function it defined p — a.e., such that 

(3.5) P(x,p)= p(x)6(p-u(x)). 

Obviously, for every fixed e > the Bohmian measure /? e is mono-kinetic by 
definition, see (|2.10p . Note however, that the limit statements for p e and J e given 
in (|3.1[) , do noi allow us to directly pass to the limit in (3 E . Thus, in general we 
can not expect the limiting Bohmian measure ft to be of the form (|3.5|) . In order 
to obtain further insight into the situation, we shall establish in the upcoming 
subsection a connection between /3 and the, by now classical, theory of Young 
measures. 
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3.2. Connection between Bohmian measures and Young measures. Con- 
sider a sequence f £ : M. d — > M m of measurable functions. Then, we recall that there 
exists a mapping p, x = p(x) : M. d — > Ai + (M. m ), called the Young measure associated 
to the sequence f £ , such that x H > (p(x),g) is measurable for all g € Co(R m ), and 
(after selection of an appropriate subsequence) : 

lim / a(x, f 6 (x)) dx = / / a(x, X)dp x (X) dx, 
£ ^° JR d Jm d jr™ 

for any function a G L^fi; C (R m )), cf. [Q Ell ES]. In view of Definition El 
we expect a close connection between the classical limit of Bohmian measures and 
Young measures. To this end, we shall first state one of the key technical lemmas 
of this work. 

Lemma 3.5. Let tp £ be uniformly bounded in H^(K. d ) with corresponding densities 
p £ , J £ G L 1 . Then, for x € M. d a.e., there exists a Young measure 

(i x : R d M + {R r x R d ), 

associated to the pair (p £ , J £ ), such that 

POO 

(3.6) (3(x,p) ^ / r d+1 p x (r,rp)dr, 



a e (x) := p £ {x)<p x, 



in the sense of measures, with equality if p £ — p in as e — > 0+. In the 

latter case, p x is a probability measure on R d+1 . 

The property of weak convergence of the particle density is crucial in order to 
express the limiting Bohmian measure j3 by (a moment of) the Young measure 
associated to p e ,J £ '. 

Proof. Assume weak convergence of p £ — p in L 1 ^), as e -> 0+. Thus, by the 
Dunford-Pettis theorem p e is uniformly integrable. Next, consider the sequence 

p £ (x) 

for if G C^°(R 2<i ;R) such that w.r.o.g. |MU~ = 1- Then \a £ \ ^ p e , and in addition 
the sequence a £ is uniformly integrable. Thus a £ — a in L l (M. ) weakly, even 
though J £ does not necessarily converge weakly in L l . In view of Definition 12 .41 we 
obviously have (l3 £ ,ip) = J a £ (x)dx, and thus also in the limit (/3,ip) = J a°(x)dx. 
On the other hand, we know that for x € M. d , the mapping 

a : (x, r, £) — > r tp ( x, — 
V r 

is continuous in r, £ and measurable in x, hence a Caratheodory function (cf. 28 ). 
From what we have seen before, we know that a{x,p £ (x),J £ {x)) = a £ (x) con- 
verges weakly in L 1 (M d ,R) and thus Theorem 2.3 in [28] asserts the existence of a 
probability measure p x , associated to (p £ , J £ ), such that 

lim J a(x, p £ (x), J £ (x)) dx = f f rip ( x,—\ dp x (r,£) dx 



*0 .lad ./si ./iJ+i \ T 

-d+1 



r p{x,p)dp x {r,rp) dx, 
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where the second line follows from the simple change of variables rp = £. Since 
p e ^ the Young measure p x has to be supported in [0, oo) x M. d and thus, we 
obtain 

= / / r d+1 tp(x,p)d^ x (r,rp)dx, 

Jm. d J[0,oo)xR d 

i.e. the assertion of the theorem, provided p s — 1 p in L 1 (R d ). 

On the other hand, if we discard the assumption of weak L 1 convergence of p e , 
we infer the existence of a Young measure p x , such that p x (R d+1 ) ^ 1, i.e. not 
necessarily a probability measure, and that (see also Proposition 4.4 in 28J): 

liminf / a(x, p £ (x), J s (x)) dx ^ / / r d+1 tp(x,p)dp x (r,rp) dx, 



E-S-0 



[0,oo)xl 



This concludes the proof. □ 

To proceed further we recall the following definition: A sequence of (measurable) 
functions {/ e }o<e^i : R d — > R is said to converge in measure to (the function) / as 
e —¥ 0, if for every 8 > 0: 

lim meas({|/ £ (x) - J{x)\ > 5}) = 0. 

Note that if in addition ^ f e e — f in Ai + (R d ) w — *, then in general: / ^ / in 
the sense of measures. 

Theorem 3.6. Letip £ be uniformly bounded in H^(R d ) with corresponding densities 

p e ,J e € L 1 . If p e — ^ p in L 1 (K d ) strongly and J E — ^ J in measure, then (3 is 
mono-kinetic, i.e. 

J(x)\ 



(3.7) P{x,p)=p{x)5[p- 



P{x) J 

and in addition J — J , where J is the measure weak — * limit established in (13.1 



Proof. We first note that strong convergence of p e in L 1 (K d ) implies that p e £ Z^f p 
in measure. Since it is known that convergence in measure of p e , J e is equivalent 
to the fact that p x is only supported in a single point of R d+1 , cf. [281 Proposition 
4.3], we conclude 

(3.8) fi x (r,0 = S(r-p(x))5(t-J(x)). 

In addition, since strong convergence of p £ in L 1 (R d ) also implies weak convergence, 
we can insert (|3.8p into (I3.6|) (with equality), to obtain 

(3(x,p) = p d+1 (x)5(p(x)p- J(x)), 

and a simple change of variable yields (|3.7[) . By computing the first moment of 
(|3.7|) w.r.t. p and keeping in mind (|3.3I) . we conclude that J = J in this case. □ 



Recalling the results of [TB], we infer that the limiting measure /3 given by (|3.7|) 
can be rewritten as 

(3(x,p) = p(x) 6 (p - u{x)) , 

where u is defined p — a.e. by (|3.2p . In this case u can be considered the classical 
limit of the Bohmian velocity field. 
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Remark 3.7. In the case where p £ — > p in measure, but not necessarily weakly 
in L 1 (K. d ), we still know that the Young measure is given by (I3.8[) . However, in 
such a situation, we can only conclude that 

/3(x,p) > p{x)6 ^p- . 

If p = 0, which can happen in principle, no information on /3 is provided. 

3.3. An alternative point of view. One might want to describe the classical limit 
of fi £ by the Young measure v x associated to (p £ ,u £ ) instead of the one associated 
to (p £ ,J £ ). However, the problem with using u x instead of p x is the fact that 
u £ := j^- is only defined p £ — a.e. and thus, we can not directly obtain a result 



for v x analogous to the one given in Theorem 13.61 Rather, we need to assume the 
existence of an appropriate extension of u £ defined on all of K d , which satisfies the 
required convergence in measure. In this case, a change of variables yields 

Thus, instead of (13. 6p we obtain 

(3.9) /3(x,p)^t / rv x (r,p)dr. 



Despite the above mentioned drawback, the measure v x is still useful to show that 
the converse statement of Theorem 13.61 is not true in general. To this end, we as- 
sume j3 to be mono-kinetic, i.e. /3(x,p) = p(x) S (p — u(x)), from which we conclude 
from $3J% that 

' p ri/x ( r, p^ — {( x >^) ^ • p = x ^ su pp p}~ 



supp : 
Thus 



supp,.^^ C {(r,p) £l M :p = u(x),r > 0} U {(r = 0,p) : p e M d } 

and we consequently infer 

Vx(r,p) = u x (r) 6(p - u(x)) + S(r)j x (p), 

where supp r oj x C (0, oo). The appearance of the second term on the right hand 
side makes the converse statement of Theorem 13.61 fail in general. In other words, 
the fact that j3 is mono-kinetic does not imply that v x is a delta distribution in p, 
which makes it impossible to conclude the strong convergence of u £ (or J £ ). This 
fact can be further illustrated by the following example. 

Example. For any ip £ £ L 2 (R d ) we can write 



(3.10) rix) = ^/pH^)e ise ^' e : 

where S £ {x) € K is defined p £ — a.e., up to additive integer multiples of 2tt. In this 
representation (which should not be confused with the WKB ansatz to be discussed 
in Section [6]) the current density reads J £ = p £ X7S £ . Assume now, that for some 
measurable set fi c K d we have 

p £ = p\ln + £>ll{Rd/n} 
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with p\ £ — in strongly and p\ £ — p 2 7^ in L 1 (IR <1 /f2) weakly. Simi- 

larly, we assume 

S £ = Sfln + SI l{m<»/n} 

with VSf £ ^°> + VSi in L°°(n) weak-* (but not strongly) and VSf VS 2 
almost everywhere on M d /f2. Then, one easily checks that 

(3(x,p) = p 2 l{ R d/o} 5{p - WS 2 (x)), 

i.e. mono-kinetic. The corresponding Young measure however, is found to be 

v x (r,p) = uj x {r)l {Rd/n} S(p - WS 2 (x)) + S(r)l n j x (p), 

where uj x (r) is the Young measure of pIlfR^/n} an d lx{p) is the Young measure 
of VSf In- In other words, the oscillations within S e do not show in the limiting 
Bohmian measure j3 (since the corresponding limiting density vanishes), but they 
do occur in the corresponding Young measure. 

Remark 3.8. It is certainly interesting to see whether the analysis given above 
directly yields information on the classical limit of the Bohmian trajectories X £ , 
P £ defined by (|2^]l . To this end, let 

T t ,x :Kt xR$^M + (M$ xlj) 

be the Young measure associated to the Bohmian trajectories (X £ (t, x), P e (t, x)) 

and assume for simplicity that p^ p Q m L 1 (M d ) strongly. Then we conclude 

from the proof of Lemma l2.5l that for all test-functions tp £ C (R X x R^), \ € Co(^t) 
it holds 

[x(t)ff ^{x 7 p)f3 £ (t,dx,dp)dt= [ X (t) f V (X e {t,x),P e (t,x))p s {x)dxdt. 
Passing to the limit e — > 0+ on both sides we find that 

(3.ii) P(t,y,p)= Y t ,x(y,p)po(x)dx. 

Formula (|3 . 1 1 1) implies that /3(t) can be uniquely determined from Tt >x but in 
general not the other way around. An immediate conclusion of (|3.11j) is that a.e. 
in t £ K : (y,p) £ supp/3(i) if and only if there exists an i £ supppo such that 
{y,p) £ suppTa-.f. In addition we infer that if 

y^ + y, r^p, 

in measure, i.e. Tt, x is only supported in a single point, then fi(t) is mono-kinetic. 
Conversely, though, from the fact that (3(t) is mono- kinetic we can only conclude 
directly that 

Tt, a = v t , x S(p- P(t,x)) 

where v± >x = ^t,x(y) is the Young measure associated to the sequence X 6 (t,x). 
We shall use (|3 . 11 1) as the basis for further investigations of the classical limit of 
Bohmian trajectories in a future work. Note, however, that the Young measure 
T t , x carries more information than is needed in order to determine the classical 
limit of Bohmian trajectories. 
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4. Comparison to Wigner measures 

In this section we shall compare the concept of Bohmian measures (and in par- 
ticular their classical limit) to the well known theory of semi-classical measures, 
also called Wigner measures, see e.g. [HI [HI HI] for a broader introduction. In the 
following, we denote the Fourier transform of a function ip(x) by 

:= f <p( x )e- ix< dx. 

4.1. Short review of Wigner measures. In order to obtain a phase space pic- 
ture of quantum mechanics one usually considers the Wigner function (or Wigner 
transformation) w £ = w £ [ip £ ], as introduced in [33] : 

(4.1) W %x,p) := £ r (, _ i y ) W ( x+ £ y ) e *P dy . 

In view of this definition, the Fourier transform of w £ w.r.t. p is given by 

(4.2) w £ (x, y) = J w(x,p)e~ iv - p dp = ^ £ (x + (a - |y) , 

and thus Plancherel's theorem together with a simple change of variables yields 

||^|| i2(R2£i) = £ ^(2^)-' i / 2 ||V/||| 2(Rt!) . 

The real- valued function w £ {t,x,p) acts as a quantum mechanical analogue for 
classical phase-space distributions. In particular, its moments satisfy 

(4.3) p £ (x) = / w £ {x,p)dp, J £ (x) = / pw e {x,p)dp, 

where the integrals on the r.h.s. have to be understood in an appropriate sense, 
since w £ i 1 (Kf X Rp) in general. 

Remark 4.1. More precisely, it is proved in [24] [19] that the Fourier transform of 
w £ w.r.t. p satisfies w £ <E Co(R^; i 1 (R^)) and likewise for the Fourier transforma- 
tion of w £ w.r.t. x 6 R d . This allows to define the integral of w £ via a limiting 
process after convolving w £ with Gaussians, cf. [24] for more details. 

The evolution equation for w £ (t, x,p) = t//[T/; e (i)] is easily derived from the linear 
Schrodinger equation (|1.1[) . It reads 

(4.4) +p- V x u> e + Q £ [V]w £ = 0, w e (0,ir,p) = w £ a (x,p), 
where Wq = w 6 ^ 6 ,} and O e [V] is a pseudo-differential operator 

(Q s [V]f)(x,p) :=-™ // ^(s, !/)/(*, (?) e^-^dydq, 
with symbol <5y £ given by 

^(^y) = i (V(W | v ) - y(x - . 

Obviously, under the assumption (IA.1I) it holds 5V £ <r Z^ ^ . V^y, in which case 
(|4.4[) formally simplifies to the classical Liouville equation on phase space. 

Note that the Wigner picture of quantum mechanics is completely equivalent to 
the Schrodinger picture. The main drawback of using w £ is that in general it can 
also take negative values and hence can not be regarded as a probability distri- 
bution. Nevertheless it has the following important property (see e.g. [19]): For 
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any operator Op 6 (a), defined by Wey I- quantization of the corresponding classical 
symbol a(x,p) £ S(M. d . x lEc), one can compute the expectation value of Op £ (a) in 
the state ip £ via 

(4.5) (ip e ,Op e (a)ip e ) L2 = // a(x,p)w e (x,p)dxdp, 

J JR 2d 

where the right hand side resembles the usual formula from classical statistical 
mechanics. To proceed further, we recall the main result proved in [24j IT9]: 

Proposition 4.2. Let ip £ be uniformly bounded in L 2 (M. d ). Then, the set of Wigner 
functions {w £ }q <£ ^i C <S'( Ri x Rp) is weak—* compact and thus, up to extraction 
of subsequences 

w £ [?P £ ] ^°> + w° = w m S'(Rt x Rp) w — *, 

where the limit w € .M + (R^ x R^) is called the Wigner measure. If, in addition 
ip £ 6 H^(M. d ) uniformly, then we also have 

p e (x) — t p(x) = / w(x,dp), J £ (x) — ? J{x) = I pw(x,dp). 

This result allows us to exchange limit and integration on the limit on the right 
hand side of (|4.5I) to obtain 

{tjf, Op e (a)-0 e )_L 2 — ? \\ a(x,p)w(x,p) dx dp. 

The Wigner transformation and its associated Wigner measure therefore are highly 
useful tools to compute the classical limit of the expectation values of physical 
observables. 

In addition it is proved in [2H [H]; that w(t) — <£> t #wo, where wo is the initial 
Wigner measure and $(£) is the classical phase space flow given by the Hamiltonian 
ODEs 

P, X(0,x,p) = x, 

■-VV(X), P(0,x,p)=p. 

In other words w(t) can be considered a weak solution of the Liouville equation. 
Note that $ t is formally obtained from (|2 . 12[) in the limit e — > 0+. It is there- 
fore natural to compare the Wigner measure associated ip £ with the corresponding 
classical limit of the Bohmian measure associated to ip £ . 

4.2. The sub-critical case. As a first step we shall prove the following basic 
result, relating j3 and w in the sub-critical case w.r.t. to the scale e. 

Theorem 4.3. Assume that ip £ is uniformly bounded in L 2 (R d ) and that in addition 

(4.7) eV^ e £ ^°> + 0, in Lf oc {R d ). 

Then, up to extraction of subsequences, it holds 

w(x,p) = f3(x,p) = p(x) S(p). 

This result can be interpreted as follows: Sequences of functions ip £ which neither 
oscillate nor concentrate on the scale e (but maybe on some larger scale), yield in 
the classical limit the same mono- kinetic Bohmian or Wigner measure with p = 0. 
Clearly, condition (|4.7[) is propagated in time by the (semi-classically scaled) free 
Schrodinger group U £ (t) = e r ltA ^ 2£ \ 



(4.6) 
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Proof. Let (p £ C^(Ri x 



ip) and write 



((3 £ , v ) - ( P E 6 p= o,v) = / p E (x)<p 



J%x) 



dx 



p £ (x) 

\7p<p (x, rj e ) ■ J £ (x) dx, 



p e (x)(p (x, 0) dx 



by using the mean value theorem. Using the fact that ip £ C^°(R^ x M~), we can 
estimate 

\(p e , ( p)-(p e 6^0,^)1 <C f \J £ (x)\ dx^C\\i/\\ L 2 {n) \\eVr\\mn), 

and since, by assumption, eVV> £ — > in L 2 oc (M. d ) we conclude that the limiting 
Bohmian measure is of the form given above. 

In order to prove the same statement for the Wigner function we again use the 
mean value theorem to write 

(4.8) r(x± S -y)=r{x)± £ -J^ Vr{x± £ -^y)-yds, 

and consider the Fourier transformation of w £ w.r.t. the variable p £ W 1 , i.e. 

w £ (x, y) = i> e (x + |y) t/j £ (x - |y) , 

as computed in (|4.2[) . Inserting (|4.8|) into w £ (x,y) we can write 

\(w £ ,p)-(p £ ,p)\ < \(R £ ^)\, 

where for every ip £ C {Ri x Ry) the remainder R £ can be estimated using the 
Cauchy-Schwarz inequality: 

l 



R e {x,y)p{x,y) dx dy 



+ yyj -v 



ds ) dx dy 



dx dy 



^e 2 c\\vr\\h 



(n)> 



where the last inequality follows from a simple change of variables. We therefore 
conclude w £ (x,y) p{x), as e — » 0+ and an inverse Fourier transformation w.r.t. 
y then yields the desired result. □ 

Remark 4.4. The proof given above, shows that the conclusion /3 — p(x) S(p) 

still holds, if (I4.7[) is replaced by the weaker assumption: J £ q, in L 1 1 oc (R d ) 

strongly. 

4.3. The case of mono-kinetic Wigner measures. In situations where we have 
concentrations or oscillation on the critical scale e the connection between /3 and w 
is much more involved. The first problem we aim to analyze in more detail, is to 
find sufficient conditions under which the limiting Wigner measure is mono-kinetic. 
We remark that mono-kinetic Wigner measures correspond to the semi-classical 
limit of quantum dynamics before caustic onset time, see |31j and Section 6 of this 
paper. 



18 



P. MARKOWICH, T. PAUL, AND C. SPARBER 



As a starting point in this direction we can state the following theorem, which can 
be seen as an analogue of Theorem l3.6l for j3. to this end, we recall the representation 
formula (|3.10[) : For any ij} £ <E L 2 ( 



4> £ {x) = ^f{x)e lSC{ - x)/e , 
with S £ {x) € K dchned p £ — a.e. (up to additive integer multiples of 2ir). 

Theorem 4.5. Let ip £ be uniformly bounded in H^(R d ), and assume p e e — °f p in 
L 1 (R ) strongly. If in addition there exists an extension of S £ to be denoted by the 
same symbol and a function S € C 1 (r2), such that 

lim ||V5 e - V5|Uoo (n)j 



+ 



where O C M. d is an open set containing supp p, then it holds 
w(x,p) = /3(x,p) = p(x) S(p — VS(x)). 

Proof. The assertion for j3 follows immediately from the dchnition of f3 £ , by using 

that J £ I p £ — \7S £ e — -4 V5 uniformly on supp p C fl- 
irt order to prove the assertion for w, we use the p-Fourier transformed Wigner 
function and the representation (|3.10p . which yields 

w s (x, y) = \Jp £ (x + y y^j \jp £ {^- y yj exp (iSS £ {x, y)) , 
where we denote the difference quotient 

(4.9) SS £ (x,y) := \ (s?{x+ \y) - S s (x- £ -y)) . 

We aim to show that w £ converges weakly to 

w(x, y) = p(x) exp (iy ■ V5(x)) . 

To this end we shall first show that ^J~pfi £ ^f m L 2 (R d ) strongly. Since, by 
assumption 

(y/ p £ (x) J dx = / p £ {x)dx £ — I p{x)dx= / (\J p(x)j dx, 

it suffices to show ^JJf — v yfp in L 2 (R d ) weakly. This, in turn, follows from a Young 
measure argument based on Proposition 4.2 and Proposition 4.3 of [28]. With the 
strong L 2 convergence at hand, we can write 

w £ (x, y) - w{x, y) = ^P £ ( x +^y) \j P £ ( x ~ \v) ~ 

+ P {X) ^( X ,y) _ el yVS( x) y 

Due to the strong convergence of y/ff and the strong continuity of the shift-operator 
on L 2 (M. d ), the first term on the right hand becomes zero in the weak limit as 
e — > 0+, i.e. after localizing with a compactly supported test-function. It remains 
to estimate the third term, for which we use 



e iSS e (x,y) _ e iyVS(x) 



< \SS £ - VS ■ y\ < \WS £ - SS\ + \SS - VS\, 
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by the mean-value theorem. Here, we first note that \SS — VS\ g, due to 

the assumed continuity of VS. For the other term we again invoke the mean- value 
theorem and write 

\VS £ - | y y V(S £ -S)(x+ y y) ds. 

It is then easily seen that the this term likewise goes to zero, as e — > 0+, by 
assumption on (the gradient of) S £ . □ 



Remark 4.6. In view of Theorem 13.61 the above given assumptions are of course 
far from optimal when one is only concerned with the limit of (3. 

Alternatively, we can also show the following variant of Theorem 14.51 where 
we impose a slightly stronger assumption on the limiting phase S. In turn, the 
assumption on the extension S £ is slightly weaker than before. 

Corollary 4.7. Let fl £ C R d be an open set containing supp p £ . Then the assertion 
of Theorem\4.5\ also holds true, if there exists an extension S £ defined on Q £ and a 



function S £ C 1 (U e ^i suc ^ ^at 



lim ||V5 e - V5|| L o. (0e) . 

e— >0 + 



Proof. The only difference from the proof given above is, that this time we write 
vT(x, y) - w(x, y) = Jp^(x + | y) Jf(x - \ y)e^ vs ^ - p(x)e^ vs ^ 



Due to the strong convergence of ^J~ff and the strong continuity of the shift-operator 
on L 2 (M. d ), the first two terms on the right hand side cancel each other in the limit 
e — » (see also Example III. 1 in [24 ). The second term can be treated similarly as 
before, using the mean- value theorem. □ 

4.4. The general case. As we have seen, we cannot expect w or /? to be mono- 
kinetic in general. It is therefore natural to study the connection between the two 
measures under more general circumstances. 

Theorem 4.8. Letip £ be uniformly bounded in H^(BL d ) with corresponding densities 
/' ■•/ • /- : . // 

£ VVF^ + 0, mLf oc (R d ) 

and if there exists an extension of S e , such that 

d 2 S £ e -kk . 

->■ 0, yi,j € l,...,d, 



e sup 



dxidxj 

where fl £ is an open set containing supp p £ , then it holds 

lim \{vf,<p)-{P s ,<p)\=0, V^Co(^xM^). 

Proof. We again consider ip £ to be given via p,10[) , and the corresponding p-Fouricr 
transformed Wigner function 



w s (x,y) = J p £ [x + ^-y) J p £ (x - ^-y ) exp (iSS e (x, y)) 
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which wc want to compare with the p-Fourier transformed representation of /3 e , 
which, in view of (|2.10[) . is given by 

fr(x,y) = p £ (x)exp(iy VS e (x)). 

To this end, we rewrite 

w £ (x,y) = ^ P £ ( x + Y y ) \f p£ { x ~Y y ) ~ P B W + P B W) exp(iSS e {x,y)), 
where SS £ is defined in (14.91) and thus 



6S e (x,y) = ± J\s'(x±^y) -yds 

= V5 £ (x) • y + j J J y T D 2 S £ (x ± ^y) y dr ds, 

with D 2 S £ denoting the Hessian matrix of S E . In other words, we have 
6S e (x,y)=VS e (x)-y + $ e (x t y) 

and thus we obtain 

ti £ (x,y) -P e {x,y) = p £ {xy WSC{x) - y (e i9 °^ - l) 



(x + ^y)Jp £ (x~^y)- P £ ( 



x) I e ivs ^>». 



In view of the assumption on S £ the first term on the right hand side goes to zero, 
as e — > 0+ and we therefore only need to take care of the second term. Using again 
the mean- value theorem we can rewrite 



+ §/ \/f(x + z)y • V 2 Up £ {x-z)) ds + p £ (x). 

I ,/n \ / z=ey/2 



ds 

z=ey/2 



Now, let ip £ Co(M 2£i ) with |supp <p\ ^ R < oo. Then, we can estimate, using the 



Cauchy-Schwarz inequality several times 

ds dx dy 



( p(x 1 y)y / p e {xT z) y ■ V z (y / p £ (x±z)) 
^C(ip)( [ p e (x)dx) ([ \eV^p~^)\ 2 dx 

\JM d J \J\x\<2R 



'\x\<2R 

where C(ip) > depends on the supp ip. By assumption, this bound goes to zero 
as e — > + , which yields the assertion of the theorem. □ 



In Subsection 15.11 we shall show that an e-oscillatory velocity field \7S £ may 
cause the limiting Bohmian measure to be different from the Wigner measure. In 
view of (|2.1ip . we also note that the condition 

eVv^^O in Lf oc (R d ), 

implies that the one part of the quantum mechanical kinetic energy which is not 
captured by the second moment of /3 £ has to converge to zero, at least locally in x. 
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In fact, it is shown in the following corollary, that this is "almost necessary" (i.e. 
at least for wave functions which are slightly more regular) to infer /3 = w. 

Corollary 4.9. Let ip £ E H^(R d ) uniformly bounded as e — > 0+ and let eVV> £ be 
compact at infinity. Furthermore assume that there exists a k > 0, such that 

H\ K+1 ^Fil) G L 2 (R d ) 

^ m L 2 (R d ). 



(4.10) 

and assume w — j3. Then eV^/Jf 
Proof. Recall that 

// \p\ 2 w £ (x, p) dxdp = 
(4.11) JjR2d 



uniformly, as e — > 0-| 



u e \x)\ 2 p £ (x) dx + e 
\p\ 2 (3 £ (x,p) dxdp - 



\W^/p £ (x)\ 2 dx 



\V^{x)\ 2 dx, 



in view of (pTTUl) and (l2"TTj) . Thus 
// \p\ 2 f3{x,p) dxdp < lim // \p\ 2 (3 £ (x,p) dx dp 



< lim 

£->0 + 



I.-'' 



\p\ w £ {x,p) dxdp 



lim 

E->0+ 



\p\ n £ (p)dp, 



where n £ denotes the momentum density, i.e. 



n £ (p) :-- 



w £ (x,p) dx 



Now, using the results given in [24j Proposition 1.7], it is easy to show that 



lim 

e->0 + 



\p\ w £ {x,p) dxdp = 



\p\ 2 w(x,p) dx dp, 



provided that tp £ satisfies the assumptions stated above. Since f3 = w, by assump- 
tion, we obtain 



lim 

E->0+ 



\p\ fi £ (x,p) dxdp 



\p\ 2 w(x,p) dx dp 



\p\ 2 (3(x,p) dx dp. 



We therefore conclude from (|4. 1 If) that eV 



in L 2 (R d ). 



□ 



In combination with Theorem 14.81 we conclude that for wave functions ip £ which 
are uniformly bounded in any e-scaled Sobolev space of higher order than H^(R d ) 

and for which e\/ip £ is compact at infinity, the fact that e\/^/ff e Zi^ o in L 2 oc (R d ), 
is indeed a necessity to obtain w = j3. 

Remark 4.10. Note that condition (|4.10j) is trivially propagated by the free 
Schrodinger dynamics corresponding to V(x) = 0. Moreover, if V(x) satisfies 
d a V € L°°(R d ) for all |a| ^ 2, a simple Grownwall estimate, combined with en- 
ergy and mass conservation, shows that (|4.10l) with k <G [0, 1] is propagated by the 
Schrodinger dynamics U £ (t) = e~ ltHe l £ on bounded time-intervals. 

5. Case studies 

In this section we shall study the case of oscillations and concentration effects on 
the (critical) scale e and compare the corresponding Bohmian and Wigner measure. 
In general we can expect all of these effects to be physically relevant, see e.g. the 
examples given in [THl 131] • 
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5.1. Oscillatory functions. Let tfj £ (x) = f(x)g(^), where / G C^(R d ;C) and 
g G C°°(R d , C) is assumed to be periodic w.r.t. some lattice L ~ Z rf , i.e. g(y + £) = 
g(y) for any y G M. d and £ 6 L. In other words, ■0 e is a slowly modulated high- 
frequency oscillation. Computing the corresponding Bohmian measure we find 

«„)-|/WP|.(f)f*(,- to (^ + .^)). 

Taking </j € Co(R^ xM^) we conclude by invoking the theory of two-scale convergence 
(see e.g. [2]), that 

^J^J Y \f(x)\"\g(y)\^^^^^ dydx, 

where Y C L is the fundamental domain of the lattice L. We thus find, that the 
limiting Bohmian measure is given by 

(5.1) P(x,p) = |/(x)| 2 ±- jf | 9 (y)| 2 5 (p - Im (^)) % 
On the other hand, let 

<?(») = E foe-**, 
e*eL* 

be the Fourier series of g(y), where L* denotes the corresponding dual lattice, and 
consider the Wigner function of %p e , after Fourier transformation w.r.t. the variable 
p G K d , i.e. 

wHx,y) = f(x+ £ -y)j(x- £ -y) E ^.e^W^^-W^)-*)). 

Then it is easy to see that, as e — >• 0+: 

^(x,y)-|/(x)| 2 E \fo\ 2 e iy - r - 

More precisely we find that 

and hence the Wigner measure associated to the ^-oscillatory function -0 e is given 

by 

(5.2) w(x,p) = \f(x)\ 2 E \h'\ 2 6(p-n 
which should be compared to (|5.1[) . 

Lemma 5.1. For g G C°°(lR d ) the limiting measures (|5.1I) and (|5.2p coincide, if 
and only if g{x) carries only a single oscillation I* G L* . 

Proof. In order to show that /3 7^ w it is enough to prove that their respective 
second moments do not coincide. To this end, we compute 

/ \p\ 2 F(x,dp) = \Y\- 1 \f(x)\ 2 [ \g(y)\ 2 

JR d JY 



dy 
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and 

/ \p\ 2 w%x,dp) = \f(x)\ 2 ]T \l 



' 2 ^i 2 = i^r 1 i/wi 2 



\Vg(y)\ 2 dy. 



Using the polar decomposition g(y) = r(y)e l6 ^ these integrals can be rewritten as 

2 



\g(y)\ 2 



Im 



\ 9(y) J 



dy= / \r(y)\ 2 \V0(y)\ 2 dy 



\Vg(y)\ 2 dy = / |Vr(y)| 2 + \r(y)\ 2 \V0(y)\ 2 dy. 



and 



Obviously these two integrals can only coincide, if |Vr(y)| = 0, which implies 
g(y ) = ce t6 ^ v \ with c ^ and 6(y) G K. In this case the support of the x- 
projection of (3 is the closure of the range of V0, i.e. bounded in "Mz. On the other 
hand, the support of the x-projection of w is L* . Hence, for a smooth function g 
the two supports can only be equal if 9{y) = y ■ £* for some t G L* , in which case 
w = p. □ 

Assume now that / is real- valued and let g(y) = e %e ^ v \ Then, the sequence ip £ 
is obviously uniformly bounded in H^(R d ) and the phase S e (x) = e0(x/s) is such 
that 

d 2 S £ d 2 9 /x\ „ . „ 



.d. 



dx e dxj dyidyj Ve> 
Therefore the first assumption of Theorem 14.81 is satisfied, but the second is not, 
unless 9 = 0. As stated above, (3 ^ w. 

5.2. Concentrating functions. We consider wave function ip £ which concentrate 
at a single point. To this end, let, for some x € R d , i(j s (x) = e~ d / 2 f with 
/ G qf(R d ;C). Thus \ip e (x)\ 2 <5(x-x ), as e ^ 0+. The corresponding Wigner 
measure has been already computed in [T8] [24] as 

(5.3) w(x, P ) = (27r)- d \f(p)\ 2 5(x-x ), 

where / denotes the Fourier transform of /. On the other hand, we easily compute 



and thus 



(/3 £ ^)= e 



-d 



(f 



S I p — Im 



tp x, Im 



/( ^) 



/I 



dx 



l/WlV ep + a; ,Im 



V /(P) y 

by a simple change of variables. We therefore conclude 

°* ' |/(p)|V(xo,to(^i) 



In other words, 

(5.4) = J(x-x ) 



dp, 



dp. 
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Again we see that the Wigner measure (|5.3|) and the classical limit of the Bohmian 
measure (15.41) are rather different in this case. 



Lemma 5.2. The limiting measures (|5 .4[) and (|5.3[) do not coincide, unless / = 0. 

Proof. Again, it is enough to prove that the second moments of /? and w do not 
coincide. By the same arguments as in the proof of Lemma f5.ll we conclude that the 
second moments can only coincide if f(y) — ce ie ^ v \ c 0, which is in contradiction 
to the fact that / € L 2 (R d ), unless c = 0. □ 

5.3. Examples from quantum physics. As a possible application we shall now 
consider some particular examples of quantum mechanical wave functions, which 
incorporate oscillatory and concentrating effects in their classical limit. 

Example (Semi-classical wave packets). In this example we consider so-called semi- 
classical wave packets (or coherent states) , which incorporate, both, oscillations and 
concentrations, i.e. 

for some given profile / € (M. d ;C). Similarly as before, we compute 



&8 e ,¥>> = e 



-d/2 



f 



X — Xq 



ip I x,po + \/elm 



V/ 



which in the limit e — > 0+ yields 

P(x,p) = \f(x)\ 2 dxS(x - x )S(p- pQ 




On the other hand, the Wigner measure of a coherent state is found in [53] to be 
w(x,p)= / \f(x)\ 2 dxS(x - x )6{p-po). 



Thus, j3 = w in this fact which makes coherent states particularly attrac- 

tive for the study of the classical limit of Bohmian dynamics [15] , Note that for 
po = this can be seen as a particular case of Theorem 14. 3\ since coherent states 
concentrate on the scale y/e. 

Example (Eigenfunctions) . Let us consider a Hamiltonian operator 

H* = ~A + V(x), 

with (real-valued) smooth confining potential V(x) +oo as \x\ — > oo. The 
corresponding spectrum is known to be discrete and the associated spectral problem 
reads 

with normalized eigenstates 4>n S L 2 (M. d ) and eigenvalues A„ € M. Now, let {e n }neN 
be a sequence such that s n n —^ Q and A^™ n — -T Aei Since V(x) is confining 
(and since ip^ 1 is normalized) there exists a subsequence, which we denote by the 
same symbol, such that 

UH 2 "^ P (^0, 
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weakly in measure. Since H 6 is self-adjoint the eigenfunctions i\) e n can be chosen 
real-valued and we therefore conclude j3{x,p) = p(x) 5{p). 

For the particular case of the harmonic oscillator V(x) = ^\x\ 2 the Wigner 
measure w, as computed in [53], is w(x,p) = 5(\x\ 2 + \p\ 2 — A), i.e. a uniform 
distribution on the energy sphere. Thus, w j3, unless A = and p{x) = S(x), 
which corresponds to a classical particle at rest, sitting at the minimum of V(x) = 
2 1 ^ | 2 - I n more generality, the fact that w ^ (3 can be concluded by invoking results 
from quantum ergodicity, see e.g. [23] or microlocal analysis, which shows that 

suppw C {x,p(E E 2d : i|p| 2 + V(x) = A}. 

6. Connection to WKB approximations 

WKB expansions are a standard approach in semi-classical approximation of 
quantum mechanics (see e.g. [THUS] and the references given therein). To this end 
one seeks an approximation of the exact solution ip e (t,x) to (|l.ip . in the following 
form 

(6.1) ^ kh (x)=a s (t,x)e iS ^ ) 

where S(t,x) G R is some e-independent (real- valued) phase function and a s (x) a 
slowly varying amplitude (not necessarily real- valued) , which admits an asymptotic 
expansion 

a e ~ a + eax + e 2 a 2 + . . . . 
Note that the ansatz (|6.ip specifies a certain e-oscillatory structure of i> e due to 
the fact that the phase Six) is assumed to be e— independent. In particular, it 
should not be confused with the representation (|3.10[) . Obviously, we find that the 
Bohmian measure in this case is given by 

P s bP^ h (t)] = \a £ (t,x)\ 2 5(p-VS(t,x)). 

Plugging (|6.1[) into the Schrodinger equation and assuming sufficient smooth- 
ness, one obtains in leading order the following equation for the the phase 

(6.2) dtS+ ^S\ 2 + V(x) =0 
and the leading order amplitude 

(6.3) d t a + Va- V5 + = 0, 

which is easily rewritten as a conservation law for p = a 2 , i.e. 

d t p + div(pVS) = 0. 

Equation (|6.2[) is nothing but the classical Hamilton- J acobi equation. If we set 
u = VS 1 , then we clearly obtain from (16.21) the inviscid field-driven Burgers equation 

(6.4) d t u+{u-V)u + VV(x) =0, 

which can formally be seen as the classical limit of (12. 7[) . 

The main problem of the WKB approach arises from the fact that (|6.2I) . or 
equivalently (|6.4I) . in general does not admit global smooth solutions. In general 
S(t, ■) G C°°(K d ) only for t G [0,T*), for some (typically small) finite time T* > 0, 
which marks the appearance of the the first caustic, or, equivalently, the appearance 
of the first shock in (|6.4[) , cf. [TO] OH]- Caustics reflect the fact that new er-scales are 
generated in the exact solution ip e (t, x )i which are no longer captured by the simple 
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ansatz ()6. 1|) . Nevertheless, at least locally in-time the WKB approximation yields 
a simple representation for ^ s (t, x) which can be extended to the case of nonlinear 
Schrodinger equations, see [TDJ [TT]. Its connection to Wigner measures has been 
extensively studied in [31 j - The connection to Bohmian measures is given in the 
following result. 

Proposition 6.1. Let Assumvtion \A.l\ and \A.S\ hold and let T* > be the caustic 
onset time. Assume there exist smooth solutions a, S € C°°([0,r*) x M. d ), with 
a(t,-) G L 2 (M. d ). Then, for the exact solution ip s (t) of the Schrodinger equation 
with WKB initial data, it holds 

/3(t,x,p) = w(t,x,p) = \a(t,x)\ 2 S(p- VS(t,x)), Vt e [0,T*). 

Proof. The statement for the Wigner measure has been proven in [16] and (in more 
generality also in [21]). In order to prove that f3 is mono-kinetic before caustics we 
refer to [T], where it is shown that for T E [0,T*): 

p eS ^ + \a\ 2 , J ee ^ + |a| 2 V5, 

in C([0, T];L 1 (R d )) strongly. Thus, recalling Theorem l3.6( we directly conclude the 
desired result. □ 

In other words, as long as the WKB approximation is valid (i.e. locally in-time 
before caustics) the classical limit of the Bohmian measure of the true solution 
to the Schrodinger equation is mono-kinetic and the same holds for the Wigner 
measure. For the latter it has been shown in [31] that locally away from caustics 
the Wigner measure can always be written as a sum of mono-kinetic terms. The 
proof requires the use of the Hamiltonian flow (14. 6|) governing w(t). Unfortunately, 
such a limiting phase space flow is not available for f3(t). All we can conclude from 
above is that for t £ [0,T*), the dynamics of (3(t) is governed by 

f X = P, X(0,x) = x, 
(6.5) { . 

{ P = -VV(X), P(0, x) = VS(0, x) = u{x). 

This is the characteristic flow associated to (|6.2p . Since it breaks down at caustics 
no information for t T* can be obtained by following this approach. In view of the 
examples given in Section [5] and the already known concentration and oscillation 
effects beyond caustics (see e.g. [THUSI]) we cannot expect a simple description of 
the classical limit of Bohmian trajectories in this case. A possible way to overcome 
this problem could be a Young measure analysis of the Bohmian trajectories in the 
spirit of [5] (see also Remark l3.8p . which, however, is beyond the scope of this work. 
We also note that in the case of of semi-classical wave packets treated in [15] , the 
problem of caustics does not appear. 

Remark 6.2. In order to give the reader a basic intuition on the limiting be- 
havior of the Bohmian measure after caustics, we recall that by stationary phase 
arguments (see e.g. [31j ) the wave function after caustics can be approximated by 
a superposition of WKB states. To illustrate the kind of phenomena which can 
happen in this situation, we consider here a sum of two WKB states, i.e. 

ip e {x) = ai e lSl( - x)/e + a 2 e iS2{x)/e , 

with real- valued oi,o 2 € C^°(R d ), S U S 2 G C°°(R d ), such that, for all x e R d it 
holds: VSi(x) ^ V^fx) and a 1 (x) > a 2 (x). 
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One the one hand, we infer from (31 j . that, in this case the Wigner measure is 
given by 

w(x,p) = a\(x) 5(p - VSi(x)) + al(x) 6{p - VS 2 {x)), 

i.e. the sum of two mono-kinetic measures. On the other hand, a lengthy but 
straightforward computation shows that the limiting Bohmian measure is given by 

1 [ 2 * 

(6.6) P(x,p) = — n(x,9)5(p-<S>(x,9))d9, 

In 



n(x, 9) := a\{x) + a 2 (x) + 2a\{x)a 2 {x) cost 



where 
and 

$ := , n . (af(xWSi(x) +(4(xWS 2 (x) + ai(x)a 2 (x) co$9CVSi(x) + VS2(x)Y) . 
n(x,0) \ J J 

To this end, we note that the computation of the current vector field Im ^ ^ ^ 

yields a smooth function which is periodic in 9(x) = (^(x) — S\{x))/e and thus 
admits a Fourier expansion w.r.t. 9. By standard two scale-convergence we infer 
that the limit as e — >• + is given by the zeroth order coefficient of this Fourier 
series, from which we deduce (I6.6[) . 

Finally, let us mention that multi-phase type WKB methods have been used 
recently, for the study of the "quantum hydrodynamic" regularisation of the Burgers 
equation [27] (see also [TBI [31] ). 
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